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Abstract 

We consider the 2D J\ — J2 classical XY model on a square lattice. In the frustrated 
phase corresponding to J2 > Ji/2, an Ising order parameter emerges by an "order due 
to disorder" effect. This leads to a discrete symmetry plus the 0(2) global one. We 
formulate the problem in a Coulomb gas language and show by a renormalization group 
analysis that only two phases are still possible: a locked phase at low temperature and 
a disordered one at high temperature. The transition is characterized by the loss of 
Ising and XY order at the same point. This analysis suggests that the 2D J\ — J2 XY 
model is in the same universality class than XY-Ising models. 
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The ground state of a large class of two-dimensional classical XY models have the partic- 
ularity to exhibit both continuous and discrete Z 2 degeneracy simultaneously in the ground 
state. It results the appearence of a new Ising-like order parameter, in addition to the con- 
tinuous U(l) symmetry It is the case of the fully frustrated XY models (FFXY) [l|]-|§, 
triangular lattice frustrated XY model (TFXY) ||, helical XY models [ 1Q| , Josephson junc- 



tions arrays in a transverse magnetic field at half flux per plaquet ||11|| . The interplay between 



both order parameters can lead to a new type of unusual critical behavior, caracterized by 
the loss of Ising and XY order at the same point H [12|]. Nevertheless, this scenario is still 



under controversy: some recent results of Olsson seem to indicate two very close transitions 
of Ising and Kosterlitz-Thouless type || or a single but decoupled transition. 
All this class of models with both continuous and discrete degeneracy are expected to be de- 
scribed by the XY-Ising model consisting of Ising and XY models coupled each other |T2], H3| 
defined by the following Hamiltonian 

h xyi = -EP+ Ba ^i) cos (^ - Q i) + G °i°i\ > (!) 

<i,j> 

where a = ±1 is the Ising spin, 6 the phase of a two component XY unit vector (XY spin) 
and < , > indicates nearest neighbors. The above Hamiltonian has a very rich phase 
structure with in particular, in the A = B plane, a branch corresponding to simultaneous 
loss of Ising and XY order |T2], 0]. The critical behavior of the FFXY and TFXY seem to 



correspond to special points of this branch. 

Nevertheless, there are some cases where the discrete Z 2 degeneracy can result of a 
continuous one which have been lifted by thermal or quantum fluctuations. It constitutes 
what is currently called an "order by disorder" effect |14], 0] , in so far as fluctuations reduce 



possible ground states. Despite the basic symmetries become identical, once fluctuations have 
been incorporated, the relations with XY-Ising model described by (|l|) do not seem a priori 
obvious. Such a situation can be encountered in spins systems with competing interactions. 
The simplest possible model displaying all these characteristics is the 2D J\ — J 2 XY model 
on a square lattice whose Hamiltonian reads 

H = -J 1 J2 cos(9i - 6j) + J 2 E cos (^ " °i) > ( 2 ) 

<i,j> «k,l» 

with Ji, J 2 > 0, and << <>> indicates the sum over next to nearest neighbors. When 
2 J\ > J 2 the ground state is ferromagnetic. It leads to a Kosterlitz-Thouless (KT) transition 



at the temperature Tkt = n ( Jl 2J2 ^ fT7|j. However, when 2Ji < J 2 , the ground state consists 



of two independent y/2 x a/2 sublattices with AF order. The ground state energy E = 
—2NJ2 does not depend on 0, an angle parametrizing the relative orientations of both 
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sublattices. This non-trivial degeneracy is lifted by thermal fluctuations and a collinear 
ordering (corresponding to = or n) is selectionned [15]. The angle <fi thus plays a role 
analogous to an Ising order parameter. Monte-Carlo simulations predict at low temperature 



a phase with nematic order and a disordered phase at high temperature [15, [Lq|. A first 



treatment which decouples both order parameters has been done in ]l7j suggesting that the 



expected KT transition is lost and replaced by another one where the Ising order parameter 
plays an important role. 

In this letter, we want to clarify the relations between the J\ — J 2 XY model and other models 
with f/(l) x Z 2 symmetry. We show that the Coulomb gas formulation of the frustrated phase 
coincides with the one of XY-Ising models when both sublattices are decoupled i.e. in the 
infrared regime. This result suggests that the Ji — J 2 XY model belongs to the class of 
XY-Ising models. Moreover, we give an estimate of the curve separating the locked and 
disordered phase in the (j^,T) plane. 

We consider the ground state where both sublattices have independent AF order. The first 
step, following Chandra et al. flBI , is to perform a gradient expansion of the classical energy 
(0). The problem is now translated in a new one on a (2 x 2) square lattice, but now with 
two spins 1 and 2 per vertices pointing in the same directions. The new classical action A 
reads 

A = -|f £ [(V0i) 2 + (V# 2 ) 2 + 2Acos0 {V x 9 l V x 9 2 - V^V^ 2 )] , (3) 



where we have defined A = ^ < 1 and introduced the lattice derivatives V x , V y |TTJ. The 



signature of the 0(2) degeneracy lies now in the strong anisotropy between x and y directions. 
The cos (ft labels the different possible classical ground states at T = 0. If we do the Gaussian 
integration, we recover the result of Henley []I5| , namely A ~ const — 0.32 ( ^j^ j proving 
that a collinear ordering (with cos</> = ±1) is effectively selected. Let us consider the action 
(§) with cos0 = 1 for example (the case cos0 = — 1 can be deduced by changing A — > —A). 
We want to include vortex excitations in the spin wave action ([$]). The usual strategy is 
to apply the Villain transformation on each quadratic terms [19[ i.e. we introduce on each 



lattice links a gauge field in order the action @ to be 2tt periodic. The drawbacks linked 
to the Villain approximation are the decoupling between spin waves and vortices. In this 
model, spin waves are responsible of the Ising order parameter appearence. The study of 
the models quoted above with both discrete and continuous degeneracy have proved this 
coupling to be strongly relevant [p]]-|J12|. To cure this drawback, we introduce, following 
Chandra et al. Hl8| , a quadrupole coupling term A c ~ — 0.32A 2 X) cos 2 (#i(r) — 62(f)), in the 



action (|3]). This term just corresponds to local spin waves effects. Hence, this quadrupole 
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coupling term plays the role of a symmetry breaking field which can be treated as follow: 



exp [h cos p(9 1 (r) — 6* 2 (r))] = exp 

S(f) 

where in our case y s = h/2 = 0.08A 2 and p = 2. To study the action (|3|), plus the symmetry 
breaking field, we first diagonalize the bilinear form in 9i and then obtain two decoupled 
actions Ai and A 2 , where the 27r periodicity is introduced in quadratic terms. The partition 
function thus reads 

z= Y. EE/ mve.e^ , (5) 

with, 

A 1 = -|E[(V%M - 27r<(r)) 2 + X[(V%(r) - 2 7 r/ 1 1 (r)) 2 - (V%(r) - 2vr/ 2 (r)) 2 f ) 

A 2 = [(V^ 2 (r) - 2nnt(r)) 2 - A[(V l fc(r) - 2^(r)) 2 - (V 2 # 2 (r) - 2vrZ 2 (r)) 2 ] 

+ tpJ2S(r)6 2 (r) + \ogy s S 2 (r) . (7) 

r 

We have introduced for each action @ and (0) four link variables per vertices nl, li and 
, l 2 with /x = 1 for x direction and \i = 2 for y direction. The n link variables correspond 
to AF interactions (proportional to J 2 ), the / link variables to ferromagnetic interactions 
(proportional to J\). It is equivalent to define two covariant derivatives on the original 
lattice, one for unit edges, one for diagonal edges. Moreover, it enables to keep traces of the 
original lattice structure, and especially of all topological excitations it can support. 
The action @ can be infered from (|7|) by taking p = and A — > —A. Because of the 2n 
periodicity, the four link variables in (0) are not independent. It still remains three degree 
of freedom per sites plus the variabla S(r). Following |T7|], we introduce three integer valued 
variables M x (r) = e^VX, M 2(r) = e^V^, M 3 (r) = VH\ - V 2 n\ (with the totally 
antisymmetric tensor), which represent three independent vortex variables. On the original 
lattice, Mi(r), M 2 (r) are associated respectively, to the circulation of the gauge field n 
around a \/2 x \/2 plaquet and of / around a 1 x 1 plaquet. M 3 (r) is a mixture of n and / 
gauge variables around a 1 x \/2 plaquet. We suppose for convenience A << 1, it will not 
change the results. After some tedious but standard manipulations, the action (^) can be 
written as 



ipS{f){6i{r) — 2 (r)) + \ogy s S 2 (r) 



(4) 
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A 2 



tt/3(1 + A)Mi(r)lof 



-M^r') + /3AM 2 (r)log 



+2/5AMi(r)log- 



^M 2 (r') - 2(3X(M 1 (r) + M 2 (r)) log • 



-M 2 (r') 
— M 3 (r') 



- A)Mi(r)6||r - r'||5(r') - ip\M 2 (r)Q\\r - r'||S(r') 



^(0 



l_ 

T ' 



+2ipM 3 (r)e||r-r'||5(r') + -f— S(r) log 

47rp 

+ ^ [log yi M?(r) + log y 2 M 2 2 (r) + \ogy 3 M 3 2 (r) + logy. S 2 (r 

r 

In this equation, we have defined a special norm 



with P 

anisotropy. In a similar way, 0||r 
fugacities defined initially by j/i = 



l-A 



+ 



l+A 



due to the 



— r 
exp( 



= arctanf^l^-^Y^). We have also included four 
7r 2 /3(l + A)/2), y 2 = exp(-vr 2 /?A/2), y 3 = 1, y s = 

0. 08A 2 exp(— p 2 /8j3) associated to the three vortices and to S(r). They can also be regarded 
as genuine variables. In order to obtain the formula (§), we have neglected terms with 
higher power of A. The coupling terms between S(r) and vortices variables can be obtained 
following ref. |2(J and using A << 1. 

Notice that M 3 (r) does not couple with itself, independently of the approximation A << 1. 
This means that vortices M 3 (r) are already present in the collinear ground state we consider 
initially To understand the role of vortices M 3 , one can study the action (|8]) with p = 
and A — > —A, in the region y^ « 1 in order to restrict the charges to take only the values 
and ±1. A straightforward generalization of Kosterlitz-Thouless equations WJ^ proves that 
the flow associated to vortices M 2 and M 3 is driven toward a high density regime as could 
have been thought (see ref. |HJ for details). Moreover A = ^j- runs quicky toward zero 

1. e the two sublattices decouple, as in the case of Heisenberg spins Hl8| . This justifies the 
approximation A << 1. A similar phenomena is naturally observed for the flow associated 
to the more general action How could we interprete this ? In fF7 |, it was argued that 
this high density regime for vortices M 2 and M 3 is responsible of the destabilization of the 
expected KT fixed point associated to Mi vortices. Following this line, one can go one step 
further. The KT flow indicates a proliferation of M 2 and M 3 vortices i.e. they try to fill the 
semi-classical vacuum, in other words they condense, so impose geometrical constraints on 
the ground state. Had we initially considered a ground state with independent AF order on 
each sublattices (so as a parameter in (§) ), the KT equations would tell us that vortices 
M 2 and M 3 condense and <fi = or 7r are the only configurations satisfying both condensate 
of vortices. Consequently, the collinear ground state can also be regarded as a lattice of 
vortex. This implies that the degrees of freedom associated to M 2 and M 3 (so to / M ) quickly 
freeze in the infrared limit. If the degrees of freedom associated to l^(r) are frozen, the 
remaining excitations are (see (§)) vortices excitations on both diagonal sublattices (Mi) 
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and local spin wave effects i. e. the symmetry breaking field, which remains relevant even in 
the regime A — ^ (see ref. ||21||). In this infrared regime the action A = A\ + A 2 considerably 
simplifies and an effective action A e ff can be drawn (where we have recoupled both actions 
and rescaled f3) 



A 



eff 



E 



■K(3J 2 M l (r) log 



l -M 1 (r')+TT/3'J 2 M[(r) log 



M[(r') 



-zp(Mi(r) + M[{r))Q\r - r'\S{r') + ^-xS(r) lo, 
+ ]T [logyi (Mi) 2 (r) + log (M() 2 (r) + lo g2/s S 2 (r) 



5(r') 



(9) 



where (3 = (3', and Mi and M{ are vortices excitations on both diagonal sublattices and 
yi, y[ their associated fugacities. The effective action flj|) corresponds to two coupled XY 
models. Under renormalization, the condition (3 = f3' is preserved and the coupling term h 
is strongly relevant and locks the phase difference in #i(r) = 9 2 {r) + kn with k = 0, 1 |]3|, |13 1 . 
It leads in the strong coupling limit (h » 1) to an effective XY-Ising model ([I]) with A = B 
and C depending on initial values of h and (3 JTS], |T2|. It follows from this analysis that 
the critical behavior of the J\ — J 2 XY model is described by a coupled doubled Coulomb 
gas which is in the class of XY-Ising models. Hence, this suggests only one transition (or 
two very closed transitions) with simultaneous loss of Ising and XY order in accordance 
with Monte Carlo simulations |15, 16]. It is worth while noting that the J\ — J 2 XY model 
enters in the universality class of XY-Ising model in a non-trivial way, by dynamical spin 
wave effects contrary to other models quoted in the introduction. The line separating the 
locked from the unlocked phase depend on the initial conditions for the fugacities. This line 
is caracterized approximately by yi,y[ ~ y s initially. It enables us to estimate the critical 
temperature T c /J 2 for a given value of A = Ji/2J 2 : 



^ ~ 0.5 * [log( 7 A 2 ) + yW^A 2 ) + 2tt 2 ] 

J2 



The associated phase diagram has been represented in Figure 1. 
T c / J 2 = 1-0 in accordance with the results of Henley [15| T c /J 2 



(10) 

For J 2 = Ji, we find 
= 0.97 ± 0.02 and of 

Fernandez et al. [IB], T c /J 2 = 0.9 ±0.02. Yet, no definitive conclusion can be drawn because 
of the lack of numerical result. In the limit J\j J 2 —>■ 0, the critical temperature converges 
logarithmically toward as it has already been checked for quantum Heisenberg spins Jl8| . 
The phase diagram of the XY-Ising model (0) consists of three branches which meet at 
a multicritical point. One of the branches separates the locked from the unlocked phase 
whereas the other two separate KT and Ising transitions. The models quoted in the intro- 
duction, where only one transition was observed seem to correspond to special points on this 
branch. Contrary to the fully frustrated XY model or triangular antiferromagnets, we have 
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one more parameter in the J\ — J2 XY model, therefore a line of transition in the ( J\j J2, T) 
plane. When J 2 >> J\, our Coulomb gas representation (|[) represents two XY models very 
weakly coupled (the central charge will be close to 2) and it is plausible to think that this 
limit corresponds to the tricritical point in the XY-Ising phase diagram where the critical 
line becomes first order. As we increase Ji, we will go away from the tricritical point in 
direction of the multicritical point where the three branches meet. So, we expect that the 
whole transition line in the (J2/J1, T) plane corresponds to a part of the critical line in the 
(A, C) plane. A similar situation was recently conjectured for a model close to FFXY model 

0- 

To sum up, we have shown using effective actions, that the J\ — J2 XY model on a 
square lattice is in the same universality class of the XY-Ising model. The transition seems 
to correspond to a whole part (between the branch and tricritical points) of the critical line 
in the phase diagram of the XY-Ising model. This could be used to test numerically the 
continuous variations of the central charge and critical exponents along this critical line. 
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Figure Caption 



FIG. 1 : The phase diagram associated to the frustrated phase of the 2D, J\ — J2 XY 
model. The transition between the disordered and the locked phases is characterized by a 
simultaneous loss of Ising and XY order like in XY-Ising models. Near the Lifshitz point, 
where the treatment fails, the curve has been extrapolated. 
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